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Abstract 

!> ! We obtain the spectra of excited heavy baryons containing one heavy quark 

00 

by quantizing the exactly-solved heavy meson bound states to Skyrme soliton. 
' The results are comparable to the recent experimental observations and quark 

model predictions, and are consistent with the heavy quark spin symmetry. 
. However, somewhat large dependence of the results on the heavy quark mass 

strongly calls for the incorporation of the soliton-recoil effects. 
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I. INTRODUCTION 



Up to the present, most ground state charm baryons containing one c-quark, from A+ 
to have been observed ]L). There have been much efforts to find excited charm baryons 
and recently the experimental evidences for £*(2530) @, A* (2593) and A* (2625) [g-§| are 
reported. Although their quantum numbers are not identified yet, the spin-parity of the 
S*(2530) is interpreted as j n = §+ and A* (2593) and A* (2625) decaying to A+7T+7T are 
regarded as candidates for j n = ~~ and |~ excited states, respectively, in accordance with the 
quark model predictions [0-0. The small mass splittings between S*(2530) and E c (2453) 
and between A* (2593) and A* (2625) are consistent with the heavy quark spin symmetry 
T3| . [T^| , according to which the hadrons come in degenerate doublets with total spin j ± 



j e ± \ (unless j e , the total angular momentum of the light degrees of freedom, is zero) in the 
limit of the heavy quark mass going to infinity. 

On the other hand, the excited heavy baryons have been extensively studied not only in 
various quark/bag models p|- p!2| , p!5| but also in heavy hadron chiral perturbation theory [16| 
and in the bound state approach of the Skyrme model pOf| ■ In the bound state model, 
the heavy baryons are described by bound states of heavy mesons and a soliton pl|j22]] . A 
natural explanation of low-lying A(1405) is one of the success of the bound state approach 
23fl , where this j n = |~ state is described by a loosely bound S'-wave K meson to soliton. 



The same picture was straightforwardly applied to the excited A*(|~) in Ref. ||17|| . The lack 
of the heavy quark symmetry in this first trial is later supplied by treating the heavy vector 
mesons on the same footing as the heavy pseudoscalar mesons p2| , and a generic structure 
of the heavy baryon spectrum of orbitally excited states is established 



However, these works were done under the approximation that both the soliton and the 
heavy mesons are infinitely heavy and so they sit on top of each other. It is evident that 
this approximation cannot describe well the orbitally and/or radially excited states due to 
the ignorance of any kinetic effects. In Ref. |19| , the kinetic effects for the excited states are 
estimated by approximating the static potentials for the heavy mesons to the quadratic form 
with the curvature determined at the origin. Such a harmonic oscillator approximation is 
valid only when the heavy mesons are sufficiently massive so that their motions are restricted 
to a very small range. The situation is improved in Ref. [20| by making an approximate 
Schrodinger-like equation and by incorporating the light vector mesons. In a recent paper 



24 1 , we have obtained all the energy levels of the heavy meson bound states by solving exactly 



the equations of motion from a given model Lagrangian without using any approximations. 
(See also Refs. p5| , p6| .) In this paper, we quantize those states by following the standard 
collective coordinate quantization method to investigate the excited heavy baryon spectra. 

In the next section, we briefly describe our model Lagrangian and the way of solving the 
equations of motion to obtain the bound states. Then, in Sec. Ill, we quantize the soliton- 
heavy-meson bound system based on the standard collective coordinate quantization method 
and derive the mass formula. The resulting mass spectra of A c , S c , A^, and baryons are 
presented in Sec. IV and compared with the recent experimental observations and with the 
quark model predictions. Some detailed expressions are given in Appendix. 
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II. THE MODEL 



We will work with a simple effective chiral Lagrangian of Goldstone mesons and heavy 
mesons [pTfl : 

(2.1) 



Here, is an effective chiral Lagrangian of Goldstone pions presented by an SU (2) matrix 
field U = exp(ir • 7r/f n ), which is simply taken as the Skyrme model Lagrangian, 



Cff = f -f Tr (d^d^U) + ^ Tr [/t^C/] 2 , (2.2) 



with the pion decay constant f n and the Skyrme parameter e. With the help of the Skyrme 
term, it supports a stable baryon- number- 1 soliton solution under the hedgehog configuration 

U (r) = exp[zr -rF(r)}, (2.3) 

where F{r) satisfies the boundary conditions, F(0) = ir and F(oo) = 0. 

The heavy pseudoscalar and vector mesons containing a heavy quark Q are represented 
by anti-isodoublet fields <3> and <£>*, with their masses m$ and m^,*, respectively. As an 
example, in case of Q = c, they are D and D* meson anti-doublets, 

$ = (7J°,7J+) and $* = (D*°,D* + ). (2.4) 

The chiral transformations of the fields are defined as 

$, $* -> $/ i t ) $*/ i t ) 1 ' ' 

where £ = VU = exp(ir • 7r/2f n ), L G SU(2) L , R G SU{2) R , and /i is an 577(2) matrix 
depending on L, R, and £. The field £ defines vector and axial vector fields as 

^ = |(^ + ^t), 
^ = f(^-^). 
Then, the covariant derivatives are expressed in terms of as 

ZVK = +Vj), (2.7) 

and a similar equation for $* which defines the field strength tensor of the heavy vector 
meson fields as = £> M $* - D u ^ . 

In our Lagrangian, we have a few parameters, /„., e, m$, m $ », /q, and qq, which in 
principle have to be fixed from the meson dynamics. We will use the experimental values for 
the heavy meson masses. In order for the quantized soliton to fit the nucleon and A masses 
the pion decay constant has been adjusted down to f n = 64.5 MeV (with e = 5.45). 
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However, recently it is shown that, taking into account the Casimir effect of the fluctuating 
pions around the soliton configuration [ffil , one can get reasonable nucleon and A masses 
with the empirical value of the pion decay constant (~93 MeV). As for the coupling constants 
fq and gq, there is no sufficient experimental data to fix them. What is known to us is that, 
in order for the Lagrangian to respect the heavy quark symmetry, they should be related to 
each other as 



lim / Q /2m$* = lim g Q = g, (2.8) 



and the nonrelativistic quark model prediction on the universal constant g is —0.75 
while the experimentally determined upper bound is g 2 <0.5 [j30"| - We will take fq and 



gq as free parameters with keeping the relation (|2.8j ) and the nonrelativistic prediction in 
mind. The parameter dependence of the results will be discussed in detail in Sec. IV. The 
Lagrangian fl2.1|) is the simplest version of the heavy meson effective Lagrangian and one 
may include the light vector meson degrees of freedom such as p and u p6| , |3Tf and the higher 
derivative terms to improve the model predictions. 

The equations of motion for the heavy mesons can be read off from the Lagrangian ( |2 . 1| ) 



as 



D^& + ml& = f Q A^, 

D^*^ + m|.$ w t = -f Q A v & + <7 Q £^ A M A $*t , 

with an auxiliary condition for the vector meson fields 

m%,D v $** = -D V D^ - ! Q D V {A V &) + 9Q^ Xp DMx%l), (2.10) 

which reduces to the Lorentz condition <9 M $*/ = in case of the free vector meson fields. 
To avoid any unnecessary complications associated with the anti-doublet structure of $ and 
$*, we work with & and $*^. When the static hedgehog configuration of Uq is substituted, 
the vector and axial vector fields become 

= (A , A) = (0, \[ai(r)r + a 2 {r)rT ■ r]), 

where v(r) = [sin 2 (F/2)]/r, ai(r) = (sinF)/r and a 2 (r) = F' — (sinF)/r. 

Then, the problem becomes to find the classical eigenmodes (especially the bound states) 
of the heavy mesons moving under the static potentials formed by the soliton field. The 
equations are invariant under parity operations and the "grand spin" rotation generated by 
the operator 

K = J + I = L + S + I, (2.12) 

where L, S, and I are the orbital angular momentum, spin, and isospin operator of the 
heavy mesons, respectively. This allows us to classify the eigenstates by the grand spin 
quantum numbers (k, k 3 ) and the parity ir. We will denote the set of quantum numbers by 
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{n} (= {n; k, k^, tt}, n is a quantum number to distinguish the radial excitations). For a 
given grand spin (k, k 3 ) with parity it = (— l) fc±1 / 2 , the general wave function of an energy 
eigenmode can be written as 

$t(r,t) = e + ^Vn(r)^g(r), <„(r,t) = e«***i<pl n (r)yfi>(i), 

where G = -i(fxL) and D^fe (?) is the (iso)spinor spherical harmonic obtained by combining 
the eigenstates of L and I. The wave functions should be normalized so that the eigenmodes 
carry a unit heavy flavor number (C = +1 and B = —1). The normalization condition is 
given in Appendix. Note the different sign convention of the energy in the exponent for the 
time evolution of the eigenmodes and that ^{r) [^(r)] is absent in case of k* 



2 _12 



Substituting Eq. ( [2.13|) into the equations of motion Q2.9Q and auxiliary condition Q2.10D 



one can obtain coupled differential equations for the radial functions, ip(r) and tp* a (r) (a = 
0, . . . , 3). (See Ref. fH for more details.) 

Given in Fig. 1 is a typical energy spectrum of bound heavy meson states obtained 
by solving the equations numerically with experimental heavy meson masses (mo = 1867 
MeV, m D * = 2010 MeV and m B = 5279 MeV, m B » = 5325 MeV), / w = 64.5 MeV, e = 5.45, 
and /q/2?ti$* = gq = —0.75. We present the binding energies defined as Ae = m $ — s. 
Comparing it with the energy spectrum obtained in the infinite heavy mass limit 



one 



can see that the "parity doubling" artifact is removed by the centrifugal energy contribution 
and there appear many radially excited states. As a trace of the heavy quark symmetry, 
the energy levels come in nearly degenerate doublets with grand spin k± = ki ± | (unless 
ki = 0) and parity n = (— l) ke [ji~8| , where = K — Sq with the heavy quark spin Sq. The 



energy levels are obtained in the soliton-fixed frame, which must be a crude approximation. 
The soliton-recoil effects should be incorporated in order for the bound state approach to 
work well with heavy flavors. In this work, however, we will proceed without incorporating 
the soliton-recoil effects. We will discuss some possible modifications of the results in Sec. 
IV, leaving the rigorous and detailed investigations to our future study. 



III. QUANTIZATION 

The soliton-heavy-meson bound system described so far does not carry any good quan- 
tum numbers except the grand spin, parity, and baryon number. In order to describe baryons 
with definite spin and isospin quantum numbers, we should quantize the system by going 
to the next order in 1/N C JO]]. This can be done by introducing collective variables to the 
zero modes associated with the invariance of the soliton-heavy-meson bound system under 
simultaneous isospin rotation of the soliton together with the heavy meson fields: 

£M) = C(tKo(r) C*(f), 

$(r,t) = $ bf (r,t)Ct(t), (3.1) 
S*(r,t) = $^(r,t)C7t(t), 

where £g = Uq and C{t) is an SU(2) matrix. The subscript "bf is to denote that they 
are the fields in the body-fixed (isospin co-moving) frame. (Hereafter, we will drop it to 
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shorten the notation and all the heavy meson fields appearing in equations are those in the 
body-fixed frame unless specified.) Assuming sufficiently slow collective rotation, we will 
work in the Born-Oppenheimer approximation where the bound heavy mesons remain in an 
unchanged classical eigenmode. 

Introduction of the collective variables as Eq. (|3.1| ) leads us to an additional Lagrangian 
ofO(l/JV c ), 

L_i = \lu 2 + u> ■ 0, (3.2) 
where the angular velocity u) of the collective rotation is defined by 

C^C = |r • w, (3.3) 



and X is the moment of inertia of the soliton [28]. The explicit form of is 



= fd^Ty^-^ + ^Ty^-^Ty^) 
+ {[^TyD^j - (D^Ty®*^] 



(3.4) 



where {A, B} + = AB + BA and 



Tv = Wo T Z° + = t x {r)r + t 2 (r)rr ■ r, 

T A = U^r^-i Q T^)^h{r){ 



(3.5) 



r x r 



with ti(r) = cosF, t 2 (r) = 1 — cosF, and h(r) = sinF. Note that it is nothing but the 
isospin current of the heavy mesons interacting with Goldstone bosons (modulo the sign) as 
discussed in Ref. |18[ . 



The spin operator J and isospin operator I of the system can be obtained by applying 
the Nother theorem to the invariance of the Lagrangian under the corresponding rotations: 

h = D ab (C)R b , 

(3.6) 

J = R+K bf , 

where R is the rotor spin conjugate to the collective variables, 

R = ^ = lev + 0. (3.7) 

Kbf is the grand spin operator of the heavy meson fields (in the isospin co-moving frame) 
and D ab (C)[= |Tr(r a Cr b C''")] is the adjoint representation of the collective variables. Note 
that the grand spin operator plays the role of the spin operator for the heavy mesons, that 
is, their isospin is transmuted into a part of the spin. 

The physical heavy baryon states with spin-parity j 71 " and isospin i can be obtained by 
combining the rotor spin eigenstates and the heavy meson bound states of grand spin k with 
the help of the Clebsch-Gordan coefficients: 
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I hh\f,h)) = J2(hJ3-h,k,k 3 \ j,j 3 ) | i;i 3 ,j 3 -k 3 } I n;k,k 3 ;7r), (3.8a) 

with k — \j — i\, \j — i\ + 1, . . . , i +j. Here, | i\ mi,m 2 } (mi, m 2 = —i, — i + 1, . . . ,i) denotes 
the eigenstate of the rotor-spin operator R a : 

K 2 \i; mi, m 2 } = i(i + 1) | i; m 1 ,m 2 }, 

R z | i;mi,m 2 } = m 2 | i\m 1 ,m 2 \, (3.8b) 
I z | i; m ll m 2 } =m 1 \i; m 1 ,m 2 } } 

and | n; k, k 3 ; 7r) is the single-particle Fock state of the heavy meson fields where one classical 
eigenmode of the corresponding grand spin quantum number is occupied: 

Km | n; k, k 3 ; vr) = k(k + 1) | n; k, k 3 ; vr), gc 
K h{ , z | n; k, k 3 ; vr) = k 3 \ n; k, k 3 ; vr). 

As an artifact of large iV c feature of the Skyrme model, one may have baryon states with 
higher isospin i > 2. We will restrict our considerations to the heavy baryons with % = 
(Aq) and i = 1 (Eg). To be precise, in Eq. ( |3.8j ), one may have to sum over the possible k 

as 

I hh]f,h)) = ^a k {i,j 3 -k 3 ,k,k 3 \ j,j 3 ) | i;i 3 ,j 3 -k 3 } \ n; k, k 3 ; it), (3.9) 

k,k3 

with the expansion coefficients to be determined by diagonalizing the Hamiltonian. How- 
ever, as far as the heavy baryon states are concerned, the mixing effects are rather small. 
It is shown in Ref. ]T8| that there is no mixing even when two states become degenerate 
in the mq — > oo limit. (Such a mixing effect plays the most important role in establishing 
the heavy quark symmetry in the pentaquark baryons [0].) Thus, we will involve only one 
single-particle Fock state in the combination ( |3.8j ). 

The physical baryons should be the eigenstates of the Hamiltonian and their masses come 
out as eigenvalues. The Hamiltonian can be obtained by taking the Legendre transformation 
with the collective variables and the heavy meson fields taken as dynamical degrees of 
freedom. Up to the order of 1/N C , we have 

H = H +1 + H° + H-\ (3.10a) 

where H m (m = +1,0, —1) is the Hamiltonian of 0(A c m ). The Hamiltonian of the leading 
order in N c is the soliton mass |28[ : 



H +1 = M soh (3.10b) 

and H° is the Hamiltonian of the heavy meson fields which yields the eigenenergy e n when 
acts on the single-particle Fock state | {n}): 

H° | {n}) = s n | {n}). (3.10c) 

Finally, the Hamiltonian of order of 1/N C arising from the collective rotation is in a form of 
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H- l = ^{K-®)\ (3.10d) 

We will take the 1/N C order term as a perturbation. Then, the mass of the heavy baryon 
state (|3.8| ) is obtained as 

m (hjn =M sol + e n + ^((i;f \ (R-&) 2 \ i;f)), (3.11) 

where e n is the eigenenergy of the heavy meson bound state involved in the construction of 
the state \i;j n )). If only one single-particle Fock state | n; k, k 3 ; n) is involved in Eq. Q3.8p, 
we can write the mass formula in a more convenient form as 

m {i , r) = M sol + e n + A + —[cj(j + 1) + (1 - c n )i(i + 1) - c n k(k + 1)]. (3.12) 

Here, we have used the Wigner-Eckart theorem to express the expectation value of © as 

(n; k, k' 3 ; tt | | n; k, k^, n) = —c n (n; k, k 3 ; n | Kbf | n\ k, k$] n), (3.13) 

which defines the "hyperfine splitting" constant c n . The explicit expressions for c n are given 
in Appendix. In evaluating the expectation value of © 2 , we have used the fact that © is 
the isospin operator (with opposite sign), which implies 

2 | n;k,k 3 ;7r) = § | n;k } k 3 ;7r). (3.14) 



IV. RESULTS AND DISCUSSIONS 

As for the Aq baryons that are constructed with the i = rotor spin state and one 
single-particle Fock state | n; k = j, k 3 ;ir), the mass formula can be further simplified as 

3 

m A Q (i) = M soi + e n + — = m N + m$ + Ae n , (4.1) 

where m N is the nucleon mass {m N = M so \ + 3/8X) and Ae n = e n — m$. Thus, the mass 
spectrum of Aq baryons is exactly the same as Fig. 1 with replacing the Ae = line by the 
niN + m$ threshold. However, the Aq spectrum obtained with the parameters of Fig. 1 (Set 
1) is not at the level of being compared with experiments. As can be seen in Fig. 1, the 
binding energy (~ 380 MeV) and the mass splitting (~ 200 MeV) between the first excited 
A* and the ground state are too small compared with the experimental values, ~ 520 and 
~ 310 MeV, respectively. However, we can easily improve the situation by adjusting the 
parameters within a reasonable range. Table I summarizes the parameter sets that we will 
examine and the parameter dependence of A c spectra are shown in Fig. 2. 

What we want to have is more deeply bound states with wider level splittings, which can 
be achieved if we have a deeper and narrower interacting potentials in the equations of motion 
for the heavy mesons. One way of obtaining such potentials in a given model Lagrangian is 
to take the empirical value for f n instead of /„. = 64.5 MeV. Since the soliton wave function 
F(r) is only a function of a dimensionless variable x = ef n r (in the chiral limit), the functions 
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ai(r) and 02 (r) appearing in potentials scale with the factor ef n . Furthermore, the soliton 
mass M so \ and the moment of inertia X come out in the form of ^8 



M sol = Mf n /e and 1 = 1 / (e 3 f„) , (4.2) 

with dimensionless quantities M and X that are independent of e and f n . If we are to have 
a correct A-N mass splitting (3/2X), we have to fix the value of e so that e 3 /„- does not 
change. This condition yields e = 4.82 when f n = 93 MeV, which implies that the soliton 
mass becomes so heavy as 1.4 GeV. We expect that the Casimir energy of fluctuating pions 
p9| can reduce it down to 0.87 GeV. In this work, for a comparison, we fix the nucleon mass 
to 940 MeV for all parameter sets. Compared with / 7r =64.5 MeV and e=5.45, ef m becomes 
1.3 times larger and thus the potential becomes deeper and narrower by the same factor. 
This is shown by the dashed line of Fig. 3. The change in f n alone (Set 2) helps the ground 
A c mass to come down to 2339 MeV with the A*-A c mass difference being 270 MeV. 

Another way of improving the results is simply to take a larger \qq\ value (with putting 
aside the experimental upper limit on \gq\ for a while), which makes the potential deeper. 
The dotted line in Fig. 3 is what obtained by varying gq and /g/2m$» to —0.92 while 
keeping / 7r =64.5 MeV (Set 3). Surprisingly, this nearly 20% change in coupling constants 
results in about 50% enhancement in the binding energy, while the mass splitting is not so 
much improved compared with that of Set 1. By the same way, we take the empirical value 
for fa and vary gq so that the ground A c mass becomes close to the experimental value, 
which is achieved with gq ~ —0.81 (Set 4). This parameter set yields comparable A c mass 
spectrum to the experiments, which looks quite encouraging. Furthermore, if we break the 
heavy quark spin symmetric relation, fq/'ZrriD* — Qq-, between the two coupling constants, 
we can obtain more realistic mass splitting between A*(|~) and A*(|~). As an example, we 
choose gq = —0.70 and fq/lm^* = —0.85 with f n = 93 MeV (Set 5). Unfortunately, these 
coupling constants are not close to the recent estimates of —0.2 ~ —0.5 We regard this 
fact as an indication of the important role of higher order corrections such as light vector 
mesons. 

In Fig. 4, we present our results on A c spectrum (obtained with parameter Set 5) together 
with the experimental values and the other model calculations; SM (Skyrme model with the 
heavy pseudoscalar mesons only) |jI7| , QM1 (quark model) J7|, QM2 ||, and QM3 |J. Our 
result can compete with the quark model calculations quantitatively. Especially, one can 
notice that it becomes much more improved compared with the first trial, SM |T7[|, in the 



Skyrme model. One may improve the result by adjusting all the parameters for the best fit. 
What we have done in this work is just to vary two coupling constants around the values 
given by the nonrelativistic quark model prediction and the heavy quark symmetric relation, 
i.e., fq/2m< s> * = gq = —0.75. As for the other parameters, we used the empirical value for 
f n with e being fixed by A- A" mass splitting, and experimental values for the heavy meson 
masses. 

Given in Fig. 5 are the spectrum of S c of our prediction (obtained with Set 5) and 
the other model calculations. As a reference line, the ground state of A c is adopted. Our 
result is again comparable to the others. However, the splitting between £* [S c (| + )] and 
S c appears too small compared with the experimental data. Note that the same is true for 
the quark models except QM2. It is also interesting to note that one cannot improve this 
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situation simply by making the hyperfine constant c larger in any way. By eliminating c n 
and I in the mass formula (|3.12|) , we obtain a model- independent relation 

m s * - m Sc = (m A - m N ) - |(m Ec - m A J. (4.3) 

(The same model-independent mass relation holds in the nonrelativistic quark model of De 
Rujula, Georgi, and Glashow @].) Thus, when our model successfully reproduces all the 
experimental values for m N , m A , m Ac , and ^ Sc , we get 

m s * - m Sc ~ 40 MeV (4.4) 

as our best prediction. It is only the half of the value evaluated with the recent S*(2533) 

We present the mass spectrum obtained with this parameter set in Fig. 6 with the 
other model calculations. The parameter set used for the charm baryons does not work well 
in the bottom sector; Set 5 yields the ground mass as 5492 MeV which is ~ 150 MeV 
below the experimental value. We may repeat the same process of varying the qq (with 
keeping the empirical value 93 MeV for f n ) to fit the mass of 5641 MeV. We also expect 
that the heavy quark symmetry relation (|2.8|) holds well in the bottom sector. This process 
leads us to qq = /g/2m B » = —0.65 (Set 6). The results with this parameter set are also 
given in Fig. 6. The mass splitting (~ 180 MeV) between the excited A£ and the ground 
state Ah appears much smaller than that of the A c given in Fig. 4, while the quark model 
calculations show nearly independent mass splittings whether the heavy constituent is a c- 
quark (~ 370 MeV) or a 6-quark (~ 330-390 MeV). Together with the differences in coupling 
constants fitting the charm baryons and the bottom baryons, this apparent difference in the 
mass splitting is certainly at odds with the heavy quark flavor symmetry. Such a heavy quark 
flavor symmetry is expected to be somehow broken because of the mass difference between 
the c-quark and 6-quark. However, since both are much heavier than the typical scale of 
the strong interaction (Aqcd ~ 200 MeV), the actual amount of the symmetry breaking in 
nature that occurs at the order of Aqcd/^q would not be so large. 

Such a behavior can be seen also in the spectrum given in Fig. 7. Since there is no 
experimental data for the baryons, we can only compare our results with the quark model 
predictions. One can find that the mass splitting between the ground Af, and the Sfe(| + ) is 
180 ~ 190 MeV, which is comparable to the quark model predictions. Also the small mass 
splitting (~ 10 MeV) between Sfe(| + ) and S fc (| + ) is still consistent with the quark model 

predictions. However, as in the A& spectrum, the excitation energy (~ 170 MeV) of Sb(| ) 
appears again smaller than the quark model values (~ 280 MeV). 

It may be the ignorance of the soliton-recoil effect in our work that causes the larger 
break down of the heavy quark flavor symmetry than what is actually implied in the model. 
In order to see this, let us go back to Fig. 1. We can see that the kinetic effect reduces 
the binding energy of the lowest D (B) meson bound state by 410 (240) MeV from its 



infinitely heavy mass limit Ae^ = -|gF'(0) ~ 790 MeV 03. Note that the ratio of the 



kinetic effects (410/240 ~ 1.7) and the ratio of energy splittings between the first excited 
state and the ground state (~ 300/200) are very close to the square root of the (inverse) 
mass ratio (a/2.6 ~ 1.6). One can easily understand this feature in the harmonic oscillator 
approximation. Thus, in our working frame, the fact that B mesons are 2.6 times heavier 
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than D mesons becomes directly reflected in the results. A simple way of estimating the 
soliton-recoil effect is to use the "reduced mass" of the soliton-heavy-meson system, as 
discussed in Refs. p3p4"f . With the soliton mass about 1 GeV, the reduced masses of the D 



mesons and B mesons become ~ 2/3 GeV and ~ 5/6 GeV, respectively. Then, the use of 
these small reduced masses can widen the energy splittings and their small ratio ~ 5/4 will 
not break the heavy quark flavor symmetry so seriously. (See also Fig. 4 of Ref. p4j.) On the 
other hand, it will require stronger potentials to overcome the larger kinetic energies, which 
should be supplied by including the light vector mesons and/or higher derivative terms into 
the Lagrangian |flj|. 

In summary, we have studied the heavy baryon spectrum in the bound state approach 
to the Skyrme model by using the exactly-solved heavy meson bound states of a given La- 
grangian. Our results are qualitatively and/or quantitatively comparable to the experimental 
observations and the quark model calculations in the charm/bottom sector. The nearly de- 
generate doublets in the spectrum are consistent with the heavy quark spin symmetry, and 
our work has a great improvement compared with the first trial ]17[ of this model. However, 
the absence of the soliton-recoil in our framework breaks the heavy quark flavor symmetry 
more than the model really implies. To be consistent with both the heavy quark spin and 
flavor symmetry, such a soliton-recoil effect should be incorporated into the picture. 
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APPENDIX: 

In this Appendix, we present the normalization condition of the heavy meson fields and 
the explicit form of the hyperfine constants. As discussed in Sec. II, the heavy meson fields 
are normalized to give a unit heavy flavor number. For a given grand spin k with parity 
7r = (— l) fc±1 / 2 , this condition can be written explicitly as 

1 = J drr 2 {2e n [\ V \ 2 + |^| 2 + A ± |^| 2 + A T |<^| 2 ] 

+ OS'Vl + ¥>iVo') + + ^Vo)^± + O0V3 + ^sVo) + 7t w ) x t 

+ I^PK + fl 2)/i T + O1V3 + ^JtVlK/^]}. ( A1 ) 

where the constants A±, fi±, and j± are written in terms of k as 
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A + = (k - l/2)(k + 1/2), \- = (k + 1/2) (A: + 3/2), 

= A: - 1/2, = -(^ + 3/2), (A2) 

7+ = H+/X+ = l/(k + 1/2), 7 _ = fi_/\_ = -l/(k + 1/2). 

The c-value defined in Eq. fl3.13|) can be written as 

(n; k, k 3 ; vr = (-l) fe±1 /2 | | n . fc> fcg . ^ = (_ 1 )*±V2) 

= -c„(n; k, k 3 ; vr = (-l) fc±1 /2 | Kfaf | ^ fc> ^ „ = (_ 1 )*±i/2 )> 

where the functionals c™, c™, and c™ are obtained by inserting Eqs. (|2.11|) , (|2.13|) , and ( p.5|) 
into Eq. fl3.4|). Their explicit expressions are as follows: 



^ = J drr 2 {e n {\^\ 2 + |^| 2 )(ti + t 2 ) + £„|^| 2 [(ti + t 2 )X ± + t 2 /i±] 

+ i(^ / Vt + ^Vs / )(ti + ^) 

+ f O0V3 + ^sVoM^Ti/^i + + ^t) [(*i + ^)A T + t 2 ^ T ]} 

- I^qC^Vs + ^sV-DOi + Oa)(*i + h)X T - \g Q (^ 2 k 2 + k\ip* 2 )t 3 ^ T 

- Iggi^h + k\{p* 2 )t z n± + \g Q {ip$k 2 + fc^s)^ 2 + /x±) 
-i9Q{<P?h + klip* 3 )h» ± }, (A4) 



c « = y drr 2 [e n \ip\H 2 -f T + e n \^\\ 2 {2t l + £ 2 ) 7t + e n \ip*\H 2 j±(l - A± - /i±) 
-i(^'Vt + ^Vo')(2ti+t 2 )7± 

+ KV0V2 + ^Vo) {^2(1 - A± - /x±) 7 i + (- + ^7t) (1 - 7±M±)} 

+ K^oVs + <ftiVo){^7±(l - 7±^±) + (- + ^7t) t( 2t i + * 2 )m=f - * 2 7±M T ]} 

- |/o(vVo + VoV)*37± - kgQ&iVl + <P*2<P*i)aiti7± 

- I^qOiVs + ^V*)ai*i(l + 7±) 

- 1 2 9qM? + l^| 2 (l + 7±)](«i + a 2 )(ti + * 2 ) 

- l9Q(<P2<P3 + 'PsVDfal + °2)0l + h)^ + ki + klipl)t 3 J± 

+ |sq(^2 + kl^Ml + 7±)(1 + /%) 

+ f^O^S + 4^2)^(1 + /^±)7± 

- ^g(^3 t/c 2 + ^2^3)^3(2 + /i±)7± 

+ |<?q(</>3% + kM)t 3 [l + /i±(l + 7± )]}, (A5) 
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c n g = J drr 2 {e n \ip\H 2l ± + e n \ip* 1 \ 2 (2t 1 + t 2 ) 1± + e n \<p* 2 \H 2 [l + 7t (1 - A± - fi±)] 
+ i(^V + ^ t O(2ti + t 2 )7± 

+ KV0V2 + ^VoH^fl + 7=f(1 - A ± - ^±)b± + (- + ^7t) 7±/^±} 
+ i(^Vs + ^Vo){^7±M± - (J + ^7t) [(2*i + + h( lTf i T - 1)]} 

+ |/o(vVo + <A)V)*37± - ^q(ViV2 + y?yl)ai*i(l + 7t) 
+ §0q(¥>iV5 + ^sVDMi7± + |^Ql^3l 2 ( a i + a 2)0i + *2)7± 

- i&jfaaVs + fsfDiai + aa)(*i + t 2 )(l ~ fJkp) - Igqi^ih + k\ip\)t^± 
~ \9Q^2^2 + kiip 2 )t 3 [-f± + fi T (l + 7±)] 

- \g Q {^h + ky 2 )h[ 1± (i + fi±) - 1] 

+ §#q(^2 + ^2^3)^3 [1 + 7±(2 + M±)] 

- f <7q(^3 + 4^)^(1 + 7±)M±}- (A6) 
The functionals fci, fc 2 j an d ^3 are defined by 



D x $*t = i{A . l(r) + k 2 (r) L^g + fc 3 (r) G^g}, (A7) 



which gives 



fcl = ~(fl ^ + ^7±^) A ± - <£>3 W A%> 

= <pl (J + ^7t) " (W + -^3) , (A8) 

k 3 = ^iv-f± - Lp 2 + -(pi 

To obtain those formulas, we have used the conjugate form of Eq. (|2.13| 
$(r,t) = e-VW^t )t ( f -) ) 

y kk 3 
i(T)t 



^(r.^-e-WWtf^), (A9) 
**(r,f) = e-« [rf(r) ^S f (r) r + ytf (r) L - ^ f (r) 3®^) (G - 2r) 



where all the operators act on the right-hand side. 
Then from Eq. (|A3|) , we can write c n as 



fc(Jfe + l)(2Jfe + l) 

+ q(n;k;7r\\G\\n;k;ir)}, 

where the "reduced matrix elements" are calculated as 



{c"(n; fc; 7r||r||n; fc; 7r) + c™(n; fc; 7r||L||n; fc; 7r) (A10) 
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(k; n = 


(-l)*+V2|| T || A . ;7r = 


( _ 1)fc+ l/ 2) 


(k; it = 


(-l) fc -V2|| T || A . ;7r = 


{-If- 1 ' 2 ) 


(k; 7r = 


(-l) fc + 1 /2|| L ||j fc;7r = 

V / II II" 


(_!)fc+l/2\ 

V / / 


(k;7i = 


(-l^-^HLll^TT = 


{-if- 1 / 2 ) 


(k;n = 


(_l)*+V2|| G || A .. 7r = 


(_ 1 )fe+l/2 ) 


(k;7i = 


(-l) fc -V2|| G || A .. 7r = 


{-if' 1 ' 2 ) 



'(fc + l)(2fc + l) 
k 



'{k + l){2k+l) 



k 



(All) 



>k{2k + l) 
k+1 



/ 2k + 1 

ik{k + iy 

I 2k + 1 
fc(fc + l)' 



and the others are zero. 



i± 



As a specific example, the c-values for the k n = ^f states are given below: 
drr 2 {§e 1+ [M 2 (-fi + t 2 ) + |tf| 2 (3fi + t 2 ) - 2\<p*\ 2 (t 1 + t 2 )\ 



Ci + 



+ K^o'Vi + vtVo )(3ti + * 2 ) + K^V; + fiVoMh + 1 2 ) 



- l^lV + a 2 ){h+t 2 )}, (A12) 



and 



Ci 



drr 2 {|ei-[|^| 2 (3ti + 1 2 ) + \<pl\ 2 {~h + t 2 ) - 2\cpl\ 2 {h + t 2 )} 

+ k^'Vi + ^;Vo )(-*i + *2) - K^oV; + ^sVo) (J - (*1 + h) 

< ?Q |<^| 2 (a 1 + a 2 )(* 1 + * 2 )}. (A13) 
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FIGURES 




FIG. 1. Energy levels of k n bound heavy meson states obtained with / w 



64.5 MeV, 



e = 5.45, m D = 1867 MeV, m D * = 2010 MeV, 



nit 



5279 MeV, m B * = 5325 MeV, and 



/(g/2m$* = qq = —0.75. The dash-dotted line is the binding energy obtained in the infinite 
mass limit. 
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FIG. 2. Parameter dependence of A c mass spectrum. 
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FIG. 5. Mass spectrum of E c (j' r ). Notations are the same as in Fig 
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FIG. 6. Mass spectrum of Ab(j n ). The predictions of Set 5 and Set 6 are presented with the 
results of SM 0, QM1 §, and QM3 |]. 
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FIG. 7. Mass spectrum of E;,^). Notations are the same as in Fi 
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TABLES 



TABLE I. Parameter sets. /„■ is in MeV unit and the others are dimensionless. 

U e gq f Q /2m$* 

Set 1 64.5 5.45 -0.75 -0.75 

Set 2 93.0 4.82 -0.75 -0.75 

Set 3 64.5 5.45 -0.92 -0.92 

Set 4 93.0 4.82 -0.81 -0.81 

Set 5 93.0 4.82 -0.70 -0.85 

Set 6 93.0 4.82 -0.65 -0.65 
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